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Purpose. To evaluate the noelin 2 gene as a disease-causing 
factor for open-angle glaucoma (OAG) and the interactions 
between the noelin 2 (OLFM2), optineurin (OPTN), and myo- 
cilin (MYOC) genes. 

Miii'iions. OI.FM2 was analyzed in 0 Japanese subjects in- 
cluding 215 patients with elevated intraocular pressure (IOP), 
2~~ wilh normal IOP, .58 with juvenile open-angle glaucoma, 
and 240 control subjects. Two single-nucleotide polymor- 
phisms (SNPs) in OPTN (c.4l2G^A and c.603T^A) and one 
SNP in MYOC (c.227G^A) were examined. Single genes were 
investigated by univariate analysis and the gene- gene interac- 
tions by logistic regression analysis. Associations between ge- 
notypes and clinical characteristics at the time of diagnosis 
were examined. 

Results. In OLFM2, 12 sequence variants were identified in 770 
Japanese subjects. Argl44Gln (exon 4) was identified in two 
(0.3%) of the patients and in none of the control subjects. 
Combinations of OLFM2/317A and OPTN /412A or OLFM2/ 
1281T and OPTN/412A were associated with patients with 
elevated IOP (P = 0.018 or P = 0.012, respectively). The 
combination of OZMf2/3 17G and OPTN/603A was signifi- 
cantly associated with elevated IOP (P = 0.018). No significant 
association was detected between SNPs in OLFM2 and 
in MYOC. Patients with normal IOP and with OLFM2/ 
678A+OPTN/412G or OZi=M2/1281C+OPTA74l2G had signif- 
icantly worse visual field scores (P = 0.022 or 0.030, respec- 
tively). 
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Conclusions. The Argl44Gln mutation in OLFM2 is a possible 
disease-causing mutation in Japanese patients with OAG. Com- 
mon polymorphisms in OLFM2 and OPT!\ may interactively 
contribute to the development of OAG, indicating a polygenic 
etiology. (Invest Ophthalmol Ms Sci. 2006; 47:5.568-5.575) 
DOI: 10. 1 l67/iovs.06-0196 

Glaucoma is the second leading cause of blindness world- 
wide, and is estimated to affect more than 60 million 
people. Open-angle glaucoma (OAG), the most common form, 
is present in almost 2% of the world's population older than 40 
years. 1 Glaucoma includes a group of conditions that is char- 
acterized by progressive optic neuropathy with visual field 
changes corresponding to the damage of the retinal nerve 
fibers. These changes are usually associated with an elevation 
of intraocular pressure (IOP). Elevated IOP is generally ac- 
cepted to be a major risk factor for glaucomatous changes. 2 

Genetic factors also play a major role in the etiology of 
OAG. 3 The first gene to be characterized was the trabecular 
meshwork-induciblc glucocorticoid response (TIGR) gene on 
lq. 4 The TIGR gene was mapped to the glaucoma locus GLC1A 
and is now known as myocilin 5 (MYOC; OMIM 601652; Online 
Mendelian Inheritance in Man; http://www.ncbi.nhn.nih.gov/ 
Omini/ provided in the public domain by the National Center 
for Biotechnology Information, Bethesda, MD). In the eye, 
myocilin is expressed in high amounts in the trabecular mesh- 
work, ciliary body, and iris and in considerably lower amounts 
in the retina and optic nerve head. 6 7 Over 50 different muta- 
tions associated witli the development of glaucoma have been 
identified in the myocilin gene in different ethnic groups 
worldwide. 8 " 18 Mutations in the myocilin gene have been 
identified in 3% to 4% of all OAG cases and in 36% of cases of 
juvenile-onset open-angle glaucoma (JOAG) in different popu- 
lations. 16 Mutations in MYOC lead to high IOP. It has been 
observed that among the three exons of MYOC, the majority 
(90?-4) of mutations are clustered in exon 3, an olfactomedin- 
like domain, whereas a few (10%) are located in exon 1, a 
myosin-like domain, and none are in exon 2. 14 Most of the 
mutations in MYOC arc located in the olfactomedin domain. 
Olfactomedin is a secreted polymeric glycoprotein of un- 
known function with an evolutionarily conserved C-terminal 
motif. 5 

Mukhopadhyay et al. 19 identified myocilin-related human 
proteins having a conserved olfactomedin domain using bioin- 
formatics approaches and examined the expression patterns in 
the eye. Myocilin-related proteins with homology to human 
myocilin were selected by BLASTp. 20 On the basis of homology 
to the N- and C termini of myocilin, three groups of proteins 
were selected: Olfactomedin 1 (noelin 1) and Pancortin iso- 
forms, OLFM2, and noelin 3 isoforms. The gene structures of 
noelin-1, -2, and -3 were determined by pair-wise BLAST 21 
analysis of the cDNA available for each gene with the retrieved 
genomic contig and located at 9q34.3, 19pl3 2, and lp22, 
respectively. Noelin 1 and 2 are expressed in the human 
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retina, 22 and noelin 3 is expressed in the human retina and in 
trabecular meshwork cells. 23 

Myocilin was found to have similar levels (60%-6l%) of 
homology, with gene products of the three noelin genes in the 
conserved olfacloniedin domains. The result of phylogenetic 
analysis indicated that myocilin may have evolved from noelin 
2 by gene duplication followed by exon fusion. 24 ' 25 The coding 
sequences of the neolin 2 gene (OLFM2) were divided into six 
exons spanning 82 kb of the genomic sequence. The gene is 
composed of a myosin-like domain (exons 1,2,3, and 4) and an 
olfactomedin-like domain (exons 5 and 6). Therefore, it is 
reasonable to investigate the olfactomedin domain containing 
proteins as potential candidates for glaucoma. 

In this study, we first screened the noelin 2 gene in Japanese 
patients with OAG for disease-causing mutations. OAG is a 
complex condition caused by multiple genes along with envi- 
ronmental factors that contribute to the phenotype. 26 " 30 In 
view of the polygenic nature of OAG, it is important to explore 
gene- gene interactions between known candidate genes. 
Therefore, we also sought to determine the possible interac- 
tions between the noelin 2, myocilin, and optineurin (OPTN; 
OMIM 602432) genes by studying single nucleotide polymor- 
phisms (SNPs) in a Japanese population. 

Methods 

Patients and Control Subjects 

A total of 770 blood samples were collected at nine institutions all over 
Japan. There were 215 patients with OAG who had elevated IOP, 277 
with normal IOP, 38 JOAG patients, and 240 individuals with no eye 
disease except cataract (control subjects). No subject was related to 
other subjects in this study. The age at diagnosis was more than 35 
years in patients with OAG and £35 years in the patients with JOAG. 
All patients received serial ophthalmic examinations, including IOP 
measurements by Goldmann applanation tonometry, Humphrey (30-2) 
or Goldmann perimetry, gonioscopy, and optic disc examination in- 
cluding fundus photographs. 

In all patients, glaucoma was diagnosed according to the following 
criteria: the presence of typical optic disc damage w ill] glaucomatous 
cupping (cup-to-disc ratio, >0.7) and loss of the neuroretinal rim; 
reproducible visual field defects compatible with the glaucomatous 
cupping: and open angles on gonioscopy. I'alienls with OAG with 
elevated IOP had an IOP >21 mm Hg at any time during the follow-up 
period. Patients with exfoliative glaucoma, pigmental-) glaucoma, or 
corticosteroid-induccd glaucoma were excluded. Patients with OAG 
who had normal IOP had glaucoma diagnosed when the untreated 
peak IOP was -21 mm I Ig at all examinations, including the three 
baseline measurements and those during the diurnal test (every 5 hours 
from 6 AM to 12 PM); when the peak IOP, with or without medication, 
alter diagnosis w as consistently -21 mm Hg throughout the follow-up 
period: and when there was an absence of a secondary cause for 
glaucomatous optic neuropathy, such as a previously elevated IOP 
after trauma, a period of steroid administration, or uveitis. 

Patients with more than -5.5 D of myopia were excluded. Four 
patients with an MYOC mutation (Ue.sOOAsn. Ala.s6.sThr, Thr448Pro, 
andPhe369Leu) 31 and tw o patients w ith an OPTS mutation (I lis26Asp)'' Q 
were also included. The procedures used in this research conformed to 
the tenets of the Declaration of Helsinki. Written, informed consent was 
obtained after the nature and possible- consequences of the study were 
explained. When applicable, the research was approved by the appropri- 
ate institutional Human Pxpcrinic illation Committee. 

The mean age at the time of blood sampling w as 647 ± 12.0 y ears 
(mean ± SD) in the patients with OAG with elevated IOP, 35.9 ± 10.7 
years in the patients with JOAG, and 60.3 ± 12.0 years in the patients 
with OAG with normal IOP. The control group was composed of 240 
volunteers who were older than 40 years and had IOPs <20 mm Hg, 
normal optic discs, and no family history of glaucoma. Mean age at the 



time of blood sampling was 69.7 ± 11.3 years. The control subjects 
were chosen to be significantly older than OAG patients with elevated 
IOP (P < 0.001), patients with JOAG (P < 0.001), and patients with 
OAG with normal IOP (P < 0.001) to try to reduce the likelihood that 
control subjects would eventually have glaucoma. 

The clinical characteristics recorded in patients with glaucoma 
w ere age at diagnosis, untreated maximum IOP (defined as IOP at 
diagnosis), and visual field defects at the initial examination (defined as 
visual field defects at diagnosis). The mean age at diagnosis was 58.0 ± 
12.1 years in patients w ith OAG with elcxated IOP. 26. 4 ± 6.0 years in 
patients with JOAG, and 57.4 ±11.5 years in patients with OAG with 
normal IOP. The mean IOP at diagnosis w as 26.7 ± 6.4 mm Hg in the 
patients with OAG with elevated IOP, 31.0 ± 10.7 mm Hg in patients 
with JOAG, and 16.6 ± 2.4 mm Hg in patients with OAG with normal 
IOP. The mean usual held score at diagnosis was 3.1 ± 0.9 in patients 
with OAG with elevated IOP, 3.0 ± 0.8 in patients with JOAG, and 
2.8 ± 0.7 in patients with OAG with normal IOP. A positive family 
history was recorded in 68 (31.6%) of 215 patients with OAG with 
elevated IOP. in 13 (34.2%) of 38 patients with JOAG, and in 86 
(31.0%) of 277 patients with OAG with normal IOP. There were 124 
(57.7%) men in the patients with OAG with elevated IOP, 27 (71.1%) 
in the patients with JOAG, 132 (47.7%) in the patients with OAG with 
normal IOP, and 114 (47.5%) in the control group. 

The severity of the visual field defects was scored from 1 to 5 
according to previously reported criteria." 83 " The data obtained by 
two types of perimeters were combined by using a five-point scale: 1, 
no alterations: 2. early defects: 3. moderate defects: 4. severe defects: 
and S. light perception only or no light perception. The first foul- 
groups on this severity scale followed Kozaki's classification based on 
Goldmann perimetry . ' ' ' 1 or the classification was based on the results 
of visual field perimetry (Humphrey Field Analyzer; Carl Zeiss Meditec, 
Dublin, CA). 35 Kozaki's classification is widely used in Japan. 

DHPLC Analysis 

Genomic DNA was isolated from peripheral blood lymphocytes by 
phenol-chloroform extraction. The six exonic coding regions of the 
noelin 2 gene were amplified by polymerase chain reaction (PCR) 
using the primer sets listed in Table I. Only exon 1 was screened by- 
direct sequencing. The other exons were screened by DHPLC analysis 
(Wave System; Transgenomic, Omaha, NH). A detailed description of 
the DHPLC technique has been published. 30 The temperatures for 
analysis of the DNA fragments are listed in Table 1 . During the high- 
throughput analysis by DHPLC, samples from three patients were 
pooled, and PCR was performed as previously described. 1 -' When the 
chromatographic patterns differed from control patterns in the pooled 
samples according to the high-throughput protocol, the samples were 
reanalyzed indi\ idually . PCR products thai show < d different chromato- 
graphic patterns from control were then sequenced. The exon 4 
region of OLFM2 was examined in 959 blood samples (770 blood 
samples and an additional 189 samples): 280 patients with OAG with 
elevated IOP. 341 patients with OAG with normal IOP. 38 patients 
with JOAG, and 300 control subjects. 

Direct Sequencing of DNA 

After purification of the PCR products, sequencing reactions were 
performed using dye termination chemistry (Prism BigDye Terminator, 
ver. 3.1 Cycle Sequencing Kit; Applied Biosystems [ABI], Foster City, 
CA), according to the manufacturer's protocol. The data were col- 
lected by a gene analy zer (Prism 310; ABI) and analyzed by computer 
(PRISM Sequencing Analysis Program, ver. 3.7; ABI). 

Genotyping SNPs 

Genotyping of c.315G^A (Alal05Ala), c.317G^A (Argl06Gln), 
c.431G^A (Argl44Gln), c.678G^A (Ala226Ala), and c.1281C^T 
(Arg427Arg) in the noelin 2 gene, and genotyping of c.412G^A 
(Thr34Thr) and c.603T^A (Met98Lys) in the optineurin gene were 
confirmed by- the chromatographic pattern of DHPLC and restriction 
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Exon Area 



Primer Sequences 



F GCGAGACCCTCACTGGGATT 
R GCCTGGAGAGGAGCTGGATT 
F GGTTGGGATTTGGGGAAGGA 
R CCAGACATGACTCCATTGTAGGAA 
F GAGTCAGAGGTTGGAGTCATGT 
R CCGTTGCTGCAGGTCCTCATA 



R GGGTGTGGCAGTCAGAGATCA 
F CCCAACTTGATCACAGCACTT 
R CTAGGCACCTATGGGCAGTCAA 
F CTAATGGCTGTAGCTGGTGCT 
R GTAGGGGAAGGTGTTGTTGTAA 
F CCAGAGCAACGTGGTGGTCA 
R CTCGTAACTGGACGTGTTGGT 
F CATGATCTGCGGTGTGCTCTA 
R GCAGCCCGAGCCACAGCATT 



:R Product 

Size PCR Tm DHPLC Tm 

(bp) (°C) (°C) 

355 is not performed 

344 67 62.0, 63.0, 64.0 

284 67 60.3, 62.3, 64.3 

249 65 62.7, 63.2, 63.7 

208 65 63.1, 64.1, 65.1 



63.7, 64.7, 65.7 
61.5, 62.0 



enzyme assay (Table 2). The G-to-A substitution at position 227 
(Met76Lys) in exon 1 of the myocilin gene w as confirmed by the chro- 
matographic pattern of 1)1 I PLC and assa\ ' (Imader assay; provided by 
the Research Department of !<&]) Center: BML. Sailama. Japan ). 

Analysis of Gene-Gene Interaction Using 
Common SNPs 

For the analysis of gene- gene interaction of MYOC and OLFM2, or 
OPTN and OLFM2 in patients with OAG, the genotypes were denned 
as follows: MYOC, c.227G^A (wild:G, mutant:A); OP 77V, c.4l2G^A 
(wild:G, mutant:A) and c.603T^A (wildT, mutant:A); and OLFM2, 
c.317G^A (wild:G, mutant:A), c.678G^A (wild:G, mutantA), and 
c.l281C^T (wild:C, mutantT) (Table 3). The interactions between 
SNPs were classified in four groups using a dominant model. The 
genotypes were MM, mutant homo/ygoie: MW, heterozygote; and 
WW, wild honiozygole. In the explanatory \ariahle. group 0 was the 
reference. 

Statistical Analysis 

Statistical analysis was performed (SPSS; SPSS Inc., Chicago, IL) and 
Hardy-Weinberg equilibrium for each SNP was determined by \ 2 test 
with 1 degree of freedom. After equilibrium was confirmed, the fre- 
quencies of the genotypes and alleles of patients w ith glaucoma were 
compared with those- of control subjects by \ J lest. Logistic regression 
analysis was used to calculate age-adjusted odds ratios (ORs) and 95% 
confidence intervals (CIs), to search for gene-gene interactions in 
patients with OAG with elevated or normal IOP. As the sample number 
of the JOAG group was insufficient to attain a normal distribution. 



patients with JOAG were excluded from the analysis. The Bonferroni 
method was used to adjust for multiple comparisons. 

Comparisons of the clinical characteristics (age, IOP, and visual 
field defect at the time of diagnosis) in patients with OAG between the 
two genotypes in a single gene or between the combined genotypes in 
the two genes were performed using the Mann-Whitney test or one- 
way ANOVA, respectively. To adjust for multiple comparisons when 
one-w a\ A\( )YA show ed a significant difference betw ecu groups (I' : 
0.05), the Tukey honestly significant difference (HSD) adjustment for 
multiple comparisons was used. P < 0.05 was considered significant. 

Bioinformatics Analysis 

To assess whether the normal residues at disease-causing mutation sites in 
the noelin 2 gene are conserved across species, the following mRNA 
entries were selected Iron i National Center for Biotechnology Information 
(NCBI; Bethesda, MD) web sites: Danio rerio (BC044164), Mus musculus 
(NM_173777), Rattus norvegicus (NM_001015017), Macaca fascicularis 
(AY650384), and Homo sapiens (NM_058164). The sequences were 
aligned by the multiple alignment tool ClustalW 37 (available online at 
http://www.ddbj.nig.ac.jp/ European Bioinformatics Institute, European 
Molecular lliolog) Laboratory Heidelberg. Germany). 

Results 

Univariate Analysis of Individual Polymorphisms 
in OLFM2 

I lie results of the sequence variants arc summarized in Table 4. 
A total of 770 Japanese subjects were studied for all exons, and 



Table 2. Genotyping OLFM2 and OPTN Sequel 



Genotype (size in bp) 



Wild Homozygote Mutant Homozygote 



C.315G ->A 
C.317G ->A 
C.431G ->A 
C.678G ->A 
C.1281C ->T 
C.412G ->A 
C.603T ->A 



C.315G (284) 
C.317G (284) 
c.431G(200, 144) 
C.678G (269) 
C.1281C (157,110) 
c.412G(188, 129) 
C.603T (277) 



c.315A(152, 132) 
c.317A(153, 131) 
C.431A (344) 
c.678A(185, 84) 
c.1281T(267) 
C.412A017) 
c.603A(175, 102) 



f/py CII t V 
Bss HII 
Hpy CH 4 IV 
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144 



human 

macaca 

rattus 

mus 

danio 



VLEQYKADT 
VLEQYKADT 
VLEQYKADT 
VLEQYKADT 
VLEQYKADA 



R TIVRLREEV 
R TIVRLREDV 
R TIVRLREEV 
R TIVRLREEV 
R V1ILRLREEV 



Figure 1. Sequence alignment of ninlanl amino acid in OLFM2 for 
Homo sapiens (human), Macaca fascial la lis (macaca). Rattus norve- 
t>icus (rallus). Mus muscuius (mns). and Danio ivrio (danio). Align- 
ment of amino acids was performed using the ClustalW alignment tool. 
Argi R)l i i was conserved among these species. 



with patients with OAG with normal IOP (P = 0.033 for allele 
frequency). 

No significant difference was detected between patients 
with glaucoma and control subjects with respect to genotypes 
or allele frequency of c.227G^A (Arg76Lys) in MYOC. 

Analysis of Gene-Gene Interactions in Patients 
with OAG 

For analysis of gene- gene interactions using the three com- 
mon SNPs in OLFM2 and one common SNP of c.4l2G->A 
(Thr34Thr) in OPTN, a logistic regression model was built 
(Table 5). A significant association was (bund in the explana- 
tory variable type 3 in this gene-gene interaction. OI.FM2/ 
317A and OPTN/412A and the OZFM2/1281T and OPTN/412A 
mutations were significantly associated with OAG with ele- 
vated IOP (P = 0.018, OR = 2.232, 95% CI: 1.150-4.331; P = 
0.012, OR = 4.240, 95% CI: 1.369-13.137, respectively). 

For analysis of gene-gene interactions using the three com- 
mon SNPs in OLFM2 and one common SNP of c.603T->A 
(Met98Lys) in OPTN, a logistic regression model was built 
(Table 6). A significant association was found in the explana- 



tory variable type 1 in this gene- gene interaction. OLFM2/non- 
31 "A and OPTN/605A were significantly associated with OAG 
with elevated IOP (P = 0.018, OR = 2.572, 95% CI: 1.173- 
5.642), and OLFM2/non-678A and OPTN/603A were weakly 
associated with OAG with normal IOP (P = 0.038, OR = 1.808, 
95% CI: 1.034 -3. 162). 

No significant interaction was detected between the three 
SNPs in OLFM2 and one common SNP of c.227G^A 
(Arg76Lys) in MYOC in patients with glaucoma (data not 
shown). 

Comparison of Clinical Characteristics of Patients 
with Glaucoma with the Combined Genotypes in 
the Two Genes 

Comparison analyses of three phenotypic variables in patients 
with OAG with elevated IOP or with normal IOP in association 
with three common SNPs of c.317G^A (Argl06Gln), 
c.678G->A (Ala226Ala), and c.1281C^T (Arg427Arg) in 
OI.FM2 did not show any statistical significance (data not 
shown). 

To examine the effect of combined genotypes in the two 
genes on considered phenotypic variables, three SNPs 
(Argl06Gln, Ala226Ala, and Arg427Arg) in OLPM2 were com- 
bined with two SNPs (Thr ViTIi and Mel98Lys) in OPTN or one 
SNP (Arg76Lys) in MYOC. Thus, nine combined genotypes in 
the two genes were obtained. The associations between the 
nine combined genotypes in the two genes and the pheno- 
typic variables (age, IOP, and visual field score) at the time 
of diagnosis of OAG in patients with elevated IOP or with 
normal IOP were determined by one-way ANOVA. Among the 
nine, the two combined genotypes of Ala226Ala+Thr34Thr 
(or c.678G^A+c.4l2G^A) and Arg427Arg+Thr34Thr (or 
c.l281C^T+c.4l2G^A) showed statistical significance in pa- 
tients with OAG with normal IOP (P = 0.036 and P = 0.039, 
respectively) in the visual field score. Patients with OAG with 
normal IOP who were OLFM2/678A and OPTN/412G carriers 
had worse visual field scores (/' = 0.022) than those who were 
OLFM2/67SG and OPTN/412A carriers based on the Tukey 



Table 5. Analysis of Gene-Gene Interactions Using Common Polymorphisms of OLFM2 and OP7jV(c.412G->A) 



SNP in Gene 




















Explanatory 


Regression 


Odds Ratio 






OLFM2 


OPTN 


Phenotype 


Variable 


Coefficient B 


Exp (B) 


95% CI 


P , 


c.317G^A 


c.412G^A 














(Argl06Gln) 


(Thr34Thr) 


OAG with elevated IOP 
OAG with normal IOP 


3 
2 


0.803 
-0.139 
0.056 
0.415 
-0.192 
-0.231 


2.232 


1.150-4.331 


0.018t 

0.862 
0.205 
0.421 
0.451 


c.678G^A 


c.412G^A 














(Aia226Ala) 


(Thr34Thr) 


OAG with elevated IOP 
OAG with normal IOP 




0.252 
-0.337 

0.435 

0.644 
-0.496 

0.042 






0.772 
0.411 
0.068 
0.381 
0.199 
0.855 


c.1281C^T 


c.412G^A 














(Arg427Arg) 


(ThtS iThr) 


OAG with elevated IOP 
OAG with normal IOP 




1.445 
0.269 
0.375 
0.621 
0.368 
0.160 


4.240 


1.369-13.137 


0.012t 
0.462 
0.128 
0.287 
0.272 
0.498 



* P by logistic regression analysis. 
tP < 0.05. 
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Table 6. Analysis of Gene-Gene Interactions Using Common Polymorphisms of OLFM2 and OP7N(c.603T>A) 
SNP in Gene 









Explanatory 


Regression 


Odds Ratio 




OLFM2 


OPTN 


Phenotype 


Variable 


Coefficient B 


Exp (B) 95°/o CI 


p* 


c.317G^A 


c.603T^A 












(Argl06Gln) 


(Met98Lvs) 


OAG with elevated IOP 


3 


0.035 




0.932 








2 


0.356 




0.127 








1 


0.945 


2.572 1.173-5.642 


0.018t 






OAG with normal IOP 




0.311 




0.369 










0.100 




0.654 
















c.678G->A 


c.603T^A 












(Ala226Ala) 


(Met98Lvs) 


OAG with elevated IOP 




-1.795 














-0.012 




0.976 








1 


0.588 




0.054 






OAG with normal IOP 




-1.008 




0.149 










0.089 




0.809 










0.592 


1.808 1.034-3.162 


0.038t 


c.1281C^T 


c.603T^A 












(Arg427Arg) 


(Met98Lvs) 


OAG with elevated IOP 




0.914 




0.122 










0.388 




0.271 










0.281 




0.377 






OAG with normal IOP 




0.161 




0.792 










0.547 




0.094 










0.511 




0.075 



HSD (Table 7). Patients with OAG with normal IOP who were 
OLFM2/1281C and OPTN/412G carriers had worse visual field 
scores (P = 0.030) than those who were OLFM2/1281T and 
OPTN/412A carriers based on the Tukey HSD. 



The phylogenelic analysis suggested thai myocilin may have 
evolved from OI.FM2 by gene duplication follow ed by exon 
fusion: The gene was composed of a myosin-like domain 
(exons 1, 2, 3, and 4) and an olfactomedin-like domain 
(exons 5 and 6). 24,25 We screened OLFM2 in Japanese pa- 
tients with OAG and normal subjects to identify disease- 
causing mutations. Twelve sequence variants were identi- 
fied: Two (Argl06Gln and Arg427Arg) have been reported 
and 10 were novel. The Argl44Gln found in exon 4 of the 
myosin-like domain was detected exclusively in two patients 
with OAG with elevated IOP and with normal IOP (0.3%, 
2/659), and was not detected in 300 normal subjects. Fingert 
et al. 14 set the criteria for probable disease-causing muta- 
tions in the myocilin gene: (1) altered myocilin amino acid 
sequence; (2) the presence of the mutation in one or more 



patients with glaucoma; (3) presence of the mutation iit less 
than 1% of the general population; and (4) absence of the 
mutation in normal individuals. The Argl44Gln mutation 
matched all four criteria. 

We also assessed the normal residues at disease-causing 
mutation sites in noelin 2 by ClustalW. Argl44 was highly 
conserved across species, including Homo sapiens, Macaca 
fascicularis, Rattus norvegicus, Mus musculus, and Danio 
rerio (Fig. 1). From the criteria for mutations and highly con- 
served sequences at the mutation site, we concluded that the 
Argl I Kiln change is a disease-causing mutation. Arginine is a 
basic amino acid and is positively charged. In contrast, glu- 
tamine is an acidic, polar uncharged amino acid. Four patients 
with MYOC mutations and two patients with OPTN mutation 
did not harbor the Argl44Gln mutation in OLFM2. There were 
no significant differences in allele frequency and genotype 
frequency between the patients with glaucoma and the normal 
subjects in the three common SNP groups of the noelin 2 gene: 
Argl06Gln, Ala226Ala, and Arg427Arg (Table 4). 

Interaction analyses of the noelin 2, myocilin, and op- 
tineurin genes were performed in patients with OAG by using 
a logistic regression model based on the combinations of com- 



Tahu: 7. Comparison of Clinical Characteristics of Patiet 



is of Two Genes 



Phenotype 


Phenotype Variable 


OLFM2/&7SG + OPTN/412A 


OLFM2/67SA + OPTN/412G 


P* 


OAG with normal IOP 


Age al diagnosis (y) 

IOP at diagnosis (mm Hg) 

Visual field score 


57.2 ± 11.3 (« = 74) 
16.5 ± 2.2 in = 70) 
2.7 ±0.6 (n = 74) 


55.2 ± 13.5 in = 17) 
15.9 ± 3.0 in = 14) 


NC$ 
NC^ 
0.022t 


Phenotype 


Phenotype Variable 


OLFM2/1281T + OPTN/412A 


OLFM2/1281C + OPTN/412G 


P* 


OAG with normal IOP 


Age at diagnosis (y) 

IOP at diagnosis (mm Hg) 


57.3 ± 7.6 (n =11) 

16.4 ± 2.2 in =11) 
2.3 ±0.5 in =11) 


57.9 ±11.5 in = 158) 
16.6 ± 2.6 in = 141) 
2.9 ± 0.7 in = 158) 


NC^ 
0.030f 



' I'll) multiple comparison (Ttiko USD). 
tP<0.05. 

$ NC, not calculated. One-way ANOVA did n 



t show a significant difference between groups CP > 0.05). 
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mon SNPs in the two genes. A significant association was found 
in the explanatory variable type 3 in two SNPs (Argl06Gln and 
Arg427Arg) in OLFM2 and Thr34Thr in OPTN, and in variable 
type 1 in two SNPs (Argl06Gln and Ala226Ala) in OLFM2 and 
Met98Lys in OPTN. These results suggest that common SNPs in 
OLFM2 and OPTN contribute interactively to OAG, indicating a 
polygenic etiology with different properties for Thr34Thr and 
Met98Lys in OPTN. 

In vivo experiments using rats, protein-protein interaction 
between optimedin and myocilin through the conserved olfac- 
tomedin domain has been demonstrated. 23 In the present 
study, however, no significant difference was detected in 
gene- gene interaction using three SNPs in OLFM2 and one 
SNP for c.22-C,^.\ (Arg _ 6I.ys) in MYOC. Arg _ 6Lys is the only 
common polymorphism reported in the Japanese population. 
If other common SNPs in MYOC are identified in the Japanese 
population, gene- gene interactions may be identified. 

Comparison analysis of three phenotypic variables of pa- 
tients with OAG with elevated IOP or with normal IOP in 
association with three SNPs of c.317G^A (Argl06Gln), 
c.678G^A (Ala226Ala), and c.1281C^T (Arg427Arg) in 
OLFM2 did not show any statistical significance. However, the 
patients with OAG with normal IOP who were OLFM2/678A 
and OPTN/412G carriers had significantly worse visual field 
scores (P = 0.022) than those who were OLFM2/67SG and 
OPTN/412A carriers. Furthermore, the patients with OAG with 
normal IOP who were OLFM2/V281C and orav/4l2G carriers 
had significantly worse visual field scores (P = 0.030) than 
those who were OZiM2/1281T and OPTN/412A carriers. 
These results suggest that the two genes may contribute inter- 
actively to clinical features in patients with O.U., indicating a 
polygenic etiology. We have reported possible gene-gene in- 
teractions between common SNPs in OPTN and TNFA for 
Japanese patients with OAG by genetic statistical analysis. 30 In 
Chinese patients with POAG, possible interactions were also 
reported between common SNPs in MYOC, OPTN, and 
APOE. 2a In view of the polygenic nature of OAG, it is quite 
important to explore possible gene- gene interactions between 
known candidate genes. This type of study will provide valu- 
able clues for further functional studies to decipher the com- 
plex pathogenesis of glaucoma. 

Recent investigations have disclosed that people classified 
as having ocular hypertension have thicker central corneal 
thickness (CCT) than do control subjects, 38 ""* 2 whereas those 
with normal-tension glaucoma have thinner CCT. 39 Unfortu- 
nately, the CCT was not measured in all subjects enrolled. The 
CCT must be considered when developing a treatment ap- 
proach for patients with ocular hypertension. 38 " 40 However, 
the epidemiologic study of glaucoma in a Japanese population, 
the Tajimi study, showed no significant intergroup difference 
in CCT among OAG with IOP > 21 mm Hg (523 ± 35 /xm), 
OAG with IOP <21 mm Hg (518±29 /mi), and subjects with- 
out glaucoma (520 ± 32 /xm), when measured by a noncontact 
specular-type instrument (SP-2000P; Topcon, Tokyo, Japan) 41 
In the normal Japanese population in the Tajimi study, the IOP 
measured with Goldman n applanation tonometry correlated 
positively with CCT, with corrected IOP (IOP reading 
-0.012 X [CCT (in /xm) - 520]) 42 This formula suggests that 
a difference in CCT of 50 /xm necessitates adjustment of IOP of 
0.6 mm Hg in the Japanese population. 

In our study, the mean IOP at diagnosis was 26.7 ± 6.4 mm 
Hg in the 215 patients with OAG with elevated IOP and 16.6 ± 
2.4 mm Hg in the 277 patients with OAG with normal IOP. For 
each patient with OAG, the CCT must be considered to make 
an exact clinical definition of the two types. Some of the 
patients with OAG with borderline IOP around 21 {n = 4) or 
22 mm Hg (« = 15) in the present study may change their 
clinical definition of the two types, depending on CCT. How- 



ever, considering that the number of patients with OAG with 
borderline IOP around 21 or 22 mm Hg was small and that 
adjustment for IOP with CCT (0.6 mm Hg/50 /xm) would be 
small in our Japanese population, the CCT would most likely 
not significantly influence our statistical analysis of association 
studies in Tables 5 and 6 in 770 subjects. 

In conclusion, the Argl44Gln mutation in the noelin 2 gene 
is a possible disease-causing mutation in Japanese patients with 
OAG. Common SNPs in OLFM2 and OPTN may interactively 
contribute to OAG, indicating a polygenic etiology. 
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